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SUMMARY 

OF THE ACTIVITY AND ITS RELEVANT RESULTS 

Some interesting results have been obtained both from spectral reflec- 
tance measurements of diseased and healthy lysimeters and compari- 
son between measured spectra and spectra calculated using the Suits ^ 
canopy model. It has been found experimentally that ‘^giallume” produces 
higher reflectance in the (0« 55 - 0c70)/un:i region of the spectrum. This 
was detected with the FG &G spectre radiometer in spite of the fact that 
the disease was not widely diffused (7% of clerotic leaves). 

The Suits ^ model was applied to the same lysimeters and all the features 
of the measured spectrum were well reproduced in the calculated one. 

As an extension of the investigation with the Suits "model a linear corre- 
lation was found between reflectance in UANDSAT band 5 a.nd the disease 
level [a disease diffusion of up to 50^), 

The specific experience of CAT A on captive low altitude balloons and 
that of ISP and JRC on poplars were combined to set up a methodology 
for performing radiometric measurements on trees from a low altitude 
R, S. platform. The technique was tested on the ISP-controlled experi- 
mental poplar farm at Frassineto. This combined use of captive balloon 
and mobile recording unit has proved itself to be an efficient; flexible 
and precise tool for tree -signature measurement. 

The processing of the considerable amount of data from the 1975 cam- 
paign on JRC lysimeters has been completed. The vegetation structure 
of rice was investigated and interpreted in dynamic terms as a signi- 
ficant factor governing the distribution of solar energy throughout the 
canopy and therefore conditions the final yield. The radiometric charac- 
teristics of the rice culture have been described for the various stages 
of development, in relation to the vegetation structure in an attempt to 
establish correlations between data of total biomass amd of grain yield. 

A preliminary study has been done over a test area (acreage roughly 
3 km ) to investigate the possibility of discriminating between rice va- 
rieties using airborne scanner data (11 channels, ground resolution 
3. 8 m). The first qualitative classification results obtained are encou- 
raging although the discrimination achieved is far from complete; im- 
provements have to be made to the sampling procedure and the atmos- 
pheric correction technique and more powerful classifiers will be tested. 
Preliminary results on the classification of natural beech forests in 
mountainous regions are also presented, A quantitative comparison be- 
tween results and ground truth is not yet available but it has been found 
that most of the beech areas are correctly identified and distinguished 
from chestnut woods present in the same area. 





1. INTRODUCTION 


The activity repox'ted here (May 15-Decsmber 31, 1976) is a continua- 
tion of the preceding one, carried out by the Joint Research Centre- 
Is.pra on the Italian test sites, following the LANDSAT-E launch, in 
coUaboraticn with the Biology Group of the Directorate General for 
Science and Education and the Italian Institutes and Organizations in 
the context of the LiANDSAT-2 follow-on investigation no, 28790 
(AGRESTE Project). 

All AGRESTE Delegation (PI and two technical Managers for the French 
and Italian test sites respectively) were invited by NASA to the GSFC to 
present the results obtained by the JRC and the collaborating Institutes 
(October 18, 1976). 

2. DANDSAT-2 IMAGERY FOR THE MADAGASCAR TEST SITE 

The situation, as regards the availability of EANDSAT-2 scenes obtained 
over our test site between 12-1-1975 and 31-7-1976, is the following; 

1. "With the exception of one, all the scenes that fulfilled the required 
conditions (i, e, overlapping the chosen area bounded by 17° $ latitude 
S ^ 18°, ^ 48° longitude E ^ 49°15'*; cloud cover ^ 60%) were re- 
ceived as black and white transparencies from ASCS Salt Lake City. 
They are listed in Table 1, with parameters used from the catalogues 
published mon thly by GSFC. 


Table 1 : Madagascar scenes available by NASA 


Observ. ID. 

Long.E 

Lat.S 

Cloud cover {%) 

Date acquired 

(a) 

2320.05595 

49.26 

17.17 

50 

•» 2/8/75 

(b) 

2375.06053 

47.34 

18.44 

40 

2/1/76 

(c) 

2392.05585 

49.21 

17.18 

• 50 

2/8/76 

(d| 

2392.05592 

49.00 

18.45 

50 

2/8/76 

(e) 

2465.06022 

48.04 

17.15 

40 

5/1/76 

(f) 

2465.06025 

47.43 

18.41 

40 

5/1/76 

(g) 

2500.05561 

49.08 

18.37 

20 

6/5/76 

(hi 

2501.06013 

, 48.04 

17.12 

20 

6/6/76 

(i) 

2501.06020 

47 ..42 

18.38 

20 

6/6/76 
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2. Positive prints of the scene received, enlarged to scale 1:700,000, 
were studied very carefully. The x'esult was twofold: 

2,1 The test area in which all the ground truth work has been done 

is located in the vicinity of the town Ambatondrazaka (a ATondra- 
zaka), the coordinates of which, are lat. 17°50"S, long. 48°25"E. 

It does not appear in the scenes labelled (b), (d), (f), (g), (i). 

2. 2 The test area is included in the scenes labelled (c), (e), (h), but 
it is entirely covered by clouds. 

3. The scene labelled (a) was not received at the JRC-Ispra. An inquiry 
was made on iebruary 15j 1977 to ASCS Salt Lake City, but from an 
examination of the picture-centre coordinates, it is doubtful that our 
test area is included in this scene. 

4. At this present stage therefore no usable scene is available. 

5. A request was made to EROS Data Centre on February 2, 1977, for 
a listing of all the scenes available which contain Ambatondrazaka, 
even with cloud cover of up to 100% and of whatever quality, acquired 
by LANDSAT-1 or -2 from the launch of LANDSAT-1 to July 31, 1976, 
All scenes available - except those alx'eady received - will then be 
ordered by the JRC as B&W transparencies, so that any possibility, 
even a very small one, to obtain usable data, will be taken into account. 

6. Coverage of our test area by LANDSAT-2 requested from and accep- 
ted by NASA for the periods of February and April/May, 1977 in the 
hope of obtaining usable data. 

The April/ May period seems to be more favourable from the point of 
view of cloud cover but the February period is necessary to achieve 
an accurate inventory of the rice (presence of water). 

There is not much sense, at this stage, in commenting on the work 
still to be done in the context of LANDSAT data-processing, in view 
. of the complete lack of data. 

3. ACTIVITY PERFORMED OVER THE RELEVANT PERIOD 
3. 1 Research Objectives and Task Distribution 

During the fourth research period of the AGRESTE Investigation, acti- 
vity has been directed towards the following objectives: 

a) Rice Investigation; 

Processing and interpretation of spectroradiometric data taken on 
rice lysimeters; 

“ to correlate reflectance with biomass and yield for di££ei*ent rice 
varieties and fertilization levels. 
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- to evaluate detectability by satellite of diseased rice under open 
•field conditions 

Development of some reference ground truth tliematic maps. 

b) Forest Investigation: 

Methodological preparation of a forest signature measurement tech- 
nique using R, S. captive balloons (test site no. 1), 

c) Computer-aided interpretation of aircraft and satellite data: 
Discrimination of rice vai*ieties from MSS aircraft data (test site nCo 1); 
Classification of natural beech forests (test site no. 3). 

Task distribution is summarized in Table 2. 

3. 2 Rice Investigation 

3.2. 1 Detectability Evaluation of Diseased Rice by Means J3f Spectro-^ 
radiometric Techniques 

3, 2, 1. 1 S pe c t r o r a diom et r ic Measurements on Some Rice 
Dy simeters 

Following the encouraging results obtained from the AGA Thermovision 
campaign'on some rice fields affected by a typical rice disease called 
*^gi?Jlume” (see 1st QPR, 2. 2.2.2) some measurements have been per- 
formed at the JRC-Ispra, in collaboration with IPW, using a continuous 
wavelength variation spectro radiometer EG&G, range (0.4 - l.O) um. 

The purpose of this experiment was the investigation of the spectral be- 
haviour of some controlled rice ly simeters v^^here the “giallume" was de- 
liberately inoculated. 

The disease was introduced into the ly simeters by means of a flying 
vector at the end of June* Radiometric observations were performed 
from July 7 to 15, when it was expected that the disease would diffuse 
in the lysimeter. 

Continuous spectral reflectances of different ly simeters on different 
days were obtained as well as reflectances in the four LANDS AT chan- 
nels. The disease did not attack all the infested lysimeters. Only one, 
presenting of clerotic leaves can be considered in the present inves- 
tigation to be really affected by '^giallume^'. 

In Fig. 1 the spectral reflectance of lysimeter no. 18, which was attack- 
ed by disease (continued line) and that of lysimeter no. 19 , which was 
not attacked by the disease (dashed line), are shown at the beginning 
(thin line) and at tlie end (thick line) of the measurement campaign. 

As a result of this investigation, it can be seen that the lysimeter in- 
fested by *’giallume” (no. 18) has a higher reflectance in the VIS ,(0. 55- 
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EBPRODUCIBILITY OF “mi; 


0^ 7) mxi region and a lower one in the NIR region than that of a non- 
infested lysimeter (no. 19 ). These results are also reported in Table 3 
where the integral reflectances on the four LANDSAT bands are shown 
for four lysimeters. of which only the no. 18 one was attacked by ^'gial- 
lume^h 

Since there was some doubt as to whether the lysimeter affected by onl^ 7% 
of clerotic leaves would be representative of the * 'gia Hume" phy to pat ho lo- 
gical behaviour for a more extended percentage of the disease for rice 
cultivation, it was decided to use the Suits ^ canopy model to obtain this 
information. 

3. 2, 1. 2 T he o r e t i c a 1 Investigation Using the Suits"' 

C ano py Model 

3.2. 1.2.1 ~ Introduction 

The 'Suits mode l(^ ) gives the directional spectral reflectance of a 
.vegetative canopy. It is based on both geometric and spectral charac- 
teristics of individual components of the considered canopy. 

Input data for the mathematical model are: 

1. The number of infinitely extended horizontal canopy layers, where 
components are randomly distributed and homogeneously mixed. 

2. The surface index of plants related to the horizontal and vertical 
components of leaves. 

3. The opticaL properties of different components of the considered 
plant (namely diffuse reflectance and transmittance of leaves, 
stalks and panicles) which are measured in the laboratory, 

3. 2. 1, 2. 2 - Comparison of Computer Evaluated Spectra and EG&G 
Measured Spectra 

The spectral reflectance for a healthy rice lysimeter, measured with 
the EG&G spectroradiometer, has been compared with the reflec- 
tance calculated by means of the Suits "'model starting from both the 
geometric and spectral characters of individual components of rice 
plants. The latter have been measured by means of a Cary 14 Lab. 
spectrometer and the results are reported in Table 4 and in Figs, 

2 and 3. Input data for the Suits ''model were measured in parallel 
with radiometric measuremexats. The sun zenith-angle of 30^, the 
angle of view of 0^ from the vertical and the azimuthal angle of 90^ 
are tlie same as those of the measurements performed using the 
EG&G spectroradiometer. 

The results of this comparison are shown in Fig. 4. 

^ Analogous comparison has been made fcr the lysimeter with 1 % of 
leaves which axe clerotic due to the effect of the disease. The results 



Table 3 : Reflectance in Landsat bands 4, 5, 6, 7 for iysimsters 17, 18, 19, 20 


7,11,15/7/75 

Band 4 

Band 5 

Band 6 

Band 7 

Band 7/ 
Band 5 

Lys.17 7.7.75 

0.0570 

0.0535 

0.2631 

0.3083 

5.7626 

11.7.75 

0.0574 

0.0528 

0.2839 

0.3398 

6.4317 

15.7.75 

0.0478 

0.0455 

0.3324 

0.4329 

9.5212 

Lys.18 7.7.75 

0.0628 

0.0570 

0.2500 

0.2935 

5.1507 

11.7.75 

0.0657 

0.0590 

0.2708 

0.3282 

5.5577 

15.7.75 

0.0534 

0.0518 

0.3036 

0.3756 

7.2562 

Lys.19 7.7.75 

0.0686 

0.0542 

0.2590 

0.3006 

5.5431 

11.7.75 

0.0590 

0.0510 

0.2867 

0.3404 

6.6713 

15.7.75 

0.0458 

0.0432 

0.3177 

0.4086 

9.4574 

Lys.20 7.7.75 

0.0595 

0.0524 

0.2453 

0.2905 

5.5464 

11.7.75 

0.0560 

0.0471 

0.2616 

0.3185 

6.7559 

15.7.75 

0.0441 

0.0431 

0.3365 

0.4264 

9.8996 


Table 4 : Geometric characterization of lysimeters N.18 and 13 (8/7/75) 


Lysimeter 

N. 18. 

— 






x = 

-14.00, 

-23.00 


S CLASS 

1 

0.05565 

0.05928 

0.05286 

0.04725 

S CLASS 

II 

0.00419 

0.00446 

0.00398 

0.00356 

S CLASS 

III 

0.0 

0.0 

0.00426 

0.01970 

S CLASS 

IV 

0:0 

0.0 

0.0 

0.0 

Lysimeter 

N. 19 







X- 

-17.00, 

-25.00 


S CLASS 

I 

0.07173 

0.09147 

0.06813 

0.07804 

S CLASS 

n 

0.0 

0.0 

0.0 • 

0.0 

S CLASS 

Ml 

U.O 

0.0 

0.00487 

0.02205 

S CLASS 

IV 

0.0 

0.0 

0.0 

0.0 

Note: 






X 



level of layers (cm) 


• CLASS 1 



horizontal and vertical 

LAI of healthy leaves 

CLASS II 



horizontal and vertical 

LAI of deseased leaves 

CLASS III 



hbriiiontal and vertical 

LAi of stalk 

CLASS IV 



free 
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are shown in Fig, 5, 

In both cases the agreement among the relevant curves is quite 
satisfactory, A.11 the features of the measui*ed spectra are very 
well reproduced in the calculated spectra. Therefore it has been 
possible to state that the Suits canopy model may be applied to eva*- 
luate the reflectance of a rice field with a high level of confidence. 


3c 2, 1, 2, 3 EKtx’apolation to Situations with more tlian 7% of Clerotic 
JLf eaves 

The model, after being checked, has been used to e^.^trapolate the 
reflectance behaviour of a diseased rice field over 7% of clerotic 
leaves. Reflectance was calculated as a function of **giaUume*‘ for 
7%, 15%, 30%, 60%, 100% of clerotic leaves. 

This study has demonstrated that the highest discrepancy in the shape 
of the spectral reflectance curves due to different level of disease 
. attack is found in the region (O, 54 - 0, 68) um. Since the band 5 of 
the LANDS AT Satellite is located approximately in this region, an 
attempt has been made to correlate between reflectance from the 
zenith in band 5 and the level of disease attack. The results of this 
correlation are shown in Fig, 6, where it may be seen that a linear - 
correlation exists up to approximately 50% of diseased leaves. 

Further investigation on open-field conditions are needed to check 
the detectability by LANDS AT satellites of ^’giallume” as a variation 
of reflectance in band 5, This hypothesis requires an epidemic dif- 
fusion of disease. 

This verification has not yet been possible owing to the absence of a 
diffuse disease attack on test site no, 1. 

3. 2, 2 ■^ppl^'^^tion of Suits Model to Some Typical Pvadiometric Problems 
for Rice 


A further extension of the Suits '"model application to tlie present study 
was made to study tlie directional reflectance of a rice field as a f\mc- 
tion of the sun- zenith angle 0, the view angle ^ , the azimuthal angle ^ 
and the characteristics of sunlight (direct or diffuse). The aim of this 
extension was to find the best -fitting conditions both for MSS flights 
and radiometric ground truth measurements. The optical and geome- 
tric properties of rice plants used in the model for the present study, 
are those measured on healthy rice lysimeters at the beginning of July , 
The following results are therefore only valid for the corresponding 
growth state of ricei 

- for a sun- zenith angle 0 20^ reflectance p of vegetative canopy 

is almost invariant as a function of the view angle ^ and the azimuthal 
'angle \p , 



refLectance(7o) 



Fig. 5 COMPARISON BETWEEN CALCULATED REFLECTANCE (?) AND RADIOMETRIC MEASUREMENT 
OF REFLECTANCE (1), DISEASE 7 ®/o 



Reflectance 



cv‘ Disease (®/o) 

Fig. 6 REFLECTANCE IN LANDSAT BAND 5 AS A FUNCTION OF DISEASE LEVEL 




- radiometric ground truth measurement may be performed with a 
view angle ^ 15^ in order" to maintain reflectance p invariant with 

a change in ^ during the day. 

reflectance p of vegetative canopy is not invariant with a change in 
sunlight characteristics (direct or diffuse). 

3. 2. 3 Reflectance Measurement Campaigns on Rice Lysimeters 

Processing of the numerous data of the 1975 campaign on JRC-lysimeters 
has been completed. It concerns three different sets of 5 lysimeters cul- 
tivated under the following conditions: 

1. dens'Ry 64 plants/m"^ at planting, per groups of 4 plantlets; i single 
addition of N-fei*tilizer (120 units), 

2. planting at the same global density but plant by plant; same N-fer- 
tilization, 

3. as in Z but N-fertilizer distributed Z/3 at planting and l/S at early 
ripening. 

The specific aim of the experiment was to investigate the vegetation 
structure as a significant factor governing the distribution of solar ener- 
gy through the canopy and hence conditioning the final yield, with an at- 
tempt to establish correlations with radiometric measurements., 

Vegetation structure is interpreted in dynamic terms as the result of 
the competition amongst tillers which occurs at different developmental 
stages of the rice plant and where two limiting factors of the environ- 
ment are mostly concerned: availability of nitrogen and radiant energy. 
Planting in groups of ple.ntlcta (treatment 1) causes an immediate and 
precocious competition between plantlets with a consequent detrimental 
effect on the number of tillers. At the end of tillering, an intense compe- 
tition between ^^major" and ^^minor^^ tillers takes place, which means 
that the further development of minor tillers can proceed only where 
there is a sufficient supply of N and light. The resulting vegetation struc- 
ture conditions the final yield. Therefore the supplementary availability 
of N at the beginning of the ripening stage (treatment 3) causes the reha- 
bilitation of tillers which would normally’ have disappeared, and a vege- 
tation structure more favourable to grain production. 

An accurate analysis of the agronomic data (dry weight and frequency 
distribution of tillers, tiller density, total leaf index and cumulative 
leaf index at several canopy layers) that aILO\VSclassi£ica.tion of the 
vegetation structure of the three treatments according to their produc- 
tion potential, with treatments 1 and 3 be mg the least and most favour- 
able, respectively. 


These differences in vegetation structure lead to significant differences 
in reflectance between the 3 treatments, both in the near IR (channel 6) 
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and in the visible light (channel 4), in spite of a similar general evo- 
- lution in time. This evolution is characterized by a maximum reflec- 
tance in the near IR 50-60 days after planting, close to the peak of the 
leaf index, and by a general decrease of the reflectance in the visible 
range during the formation of the vegetation structure, until a minimum 
is reached at the end of earing (about 80 days after planting). The ratio 
P6/P4 exhibits tire typical evolution described in earlier reports with a 
peak occurring in the middle of earing around 70 days after planting. 

A rough correlation can be found between these peak values of the ratio 
f (,/ P 4 and the respective grain yield of treatments 3 and 2 (542 and 
g/m^ respectively). The displacement in time of the peak value of 
the ratio in treatment 1 does not allow generalization of such a correla- 
tion (yield as low as 436 g/m^ as a result of a less favourable vegetation 
structure). On the other hand, the limited number of 5 duplicates inside 
each treatment is insufficient to establish a significant correlation with 
the grain yield. It must be pointed out that no single correlation between 
radiometric data and grain yield is likely to be valid without taking into 
account the analysis of the vegetation structure and of its development 
in relation to the limiting environmental factors. The variable availabi- 
lity of N in the lysimeter experiments illustrates this point very clearly 
if the availability of N is not excessive and occurs at an adequate stage 
of development (at the beginning of the ripening stage as in treatment 3 ), 
the potential of production v/ill be oriented towards an increased total 
biomass, positively correlated to the radiometric data and to the final 
grain yield. On the contrary, excess of N or an inadequate vegetation 
structure will cause a decrease of the yield coefficient (ratio of the dry 
weight of panicle to the dry weight of total biomass), thus orienting the 
productive potential towards production of vegetative material, causing 
at the same time a shift in the correlation of the grain yield with radio- 
metric data. 

,3.2.4 Ground Truth Preparation 

The ground truth map to the south of Mortara (see 3rd QPR, 3. 2.4) has 
been improved, in collaboration with ENR, in order to report on some 
additional interesting targets (water logged meadows, corn, poplars, 
grass, etc.) for computer-aided classification purposes. Eleven differ- 
ent rice varieties have already been distinguished. This map will be- 
come the foundation document for the discrimination studies on rice va- 
rieties, using Bendix airborne scanner data (see 3.4, 3), 

Work on drawing-up a ground truth map has reached an advanced stage, 
but only on rice as a whole, and on a much larger area between Mortara 
and the Po-river (test site no. 1), 
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3,3 Forest Investigation 

3. 3. 1 Metliodoloj_ical_Se^^^^^^ Signature Measurement 

Captive Balloons 

During the present intensive phase of AGPESTE research on forest 
classification, it has been found to be useful to have exhaustive infor- 
mation on the spectral signature characteristics of poplar clones and 
some inventory-concurrent tree species. The need arose for a low al- 
titude R.S. platform, in order to perform radiometric measurements 
on trees, more specifically on poplar, willow and robinia afforestations. 

Since low-altitude captive balloons represent an easily controlled and 

platform, it was decided to combine the siDecific experience 
of CATA and the ISP' and JRC ' investigation on poplars. 

Wife this purpose in mind ah experiment was performed in the autumn 
of 976 wife tiie aim of checldng and setting up this technique at Casale 
onferrato (test site no. 1) on the poplar experimental farm of Frassi- 
neto, controlled by ISP. 

A Delacoste type V. D. 4 captive balloon, equipped by CATA-s research- 
with in Exotech mod. 100 A radiometer and a remotely motor-con- 
trolUd Olympus camera for ground truth documentation was positioned 
at about 80. m above poplar experimental fields. 

Radiometric data of reflected irradiance have been sent by cable from 
the vertically suspended Exotech mod. 100 A to the JRC 's mobile unit 
on ground. These data were recorded in a Schlumberger lO-channel 

1 ™ parallel with the data coming from another Exotech mod. 

lUU A measuring the incoming sun irradiance. 

Spectral signature of different poplar clones for LAWDSAT bands were 
automafecally obtained by making a simple ratio of the recorded Exotech 
mo . A data. In this way the joint use of captive balloon and mobile 
recording unit has shown the entire R. S. system to be an efficient, fiex- 
1 > e and precise tool for tree signature measurement. 

Owing to both fee adverse meteorological conditions and the rather ad- 
vanced phenological state of poplars it was decided not to make a regu- 
lar measurement campaign on the ISP experimental area of Frassineto. 

A cornplete series of measurements is envisaged in spring 1977. Spec- 
tral signature will be investigated forj 

t 

a) Populus alba, 

b) Populus deltoide, 

c) Populus X euroamericana, 

d) Saiix alba. 
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e) Robinia, 


wbicli are the three wide clonal classes of poplar and the concurrent 
arboreous varieties in the AGRESTE Italian test sites. 

3, 3. 2 955??d_T5?!h.?Eeparatipn_fo_r Classification 

^;°^;^^^;®ntation has been prepared by INPL. NIR Aerial 
p (scale 1:17, 000) taken in 1974 have been interpreted. Some 
homogeneous zones of beeches, chestnuts and meadows have been 

Maritime 

Alps ^test site no. 3)> i, e. : 

a) Valle Stura, between Demonte and Aisone, 

b) Vallono dell^Arma, between Demonte and Rifugio Viridio. 

3.4 Satellite and Air craft Data Processing and Interpretation 
3.4.1 Implementati^ Software Package 

3, 1^1 General 

^ parallel with the mam classification and interpretation activity, the 
eve opmeni. o the interactive software package has proceeded during 
e period under consideration. A set of new modules has been created 
capable of performing the following tasks; ’ 

- extraction of statistical data from training areas, 

- supe.rvised clarification according to 3 different methods, 

~ ratios of two linear combinations of selected channels. 

- file management progmms for storage and retrieval of information 
created or needed during a session^ 

Most of these programs are in operative condition and some have reach- 
ed test status. Much care has been taken to obtain a good definition of 
an data involved in the complete system of interactive program.s, in 
or er to ensure flexible cross-over of information between different 
procedures, ease of extension of existing modules, logical construc- 
tion of foreseen algorithms and the possibility of the creation of other 
software interfaces for terminals still to be defined. 

3.4. 1.2 Colour CRT Display 

since July 1976 a COMTAL 8000 eerl.e Image procseing device tae been 

^ ^ 1 ^ configuration comprises the colour CRT Video, a 

refresh disk memory for three images of 5] 2 x 512 bytes and for three 
graphic overlays (1 bit for each point), trackball, control unit with three 
unction memories and a pseudo-colour memory, DEC PDP ll/05 pro- 
cessor with 16 Kbytes of core memory, 45 ips/800 bpi 9 -track tape unit. 
Since the software delivered with the COMTAE unit is rather limited and 
in particular lacks a true operating system, some additional routines 
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in particular one for tape input written in the scanner data format 
used at the JRC, had to be added in machine code of the PDP 11* Also 
in die batch program MAPEDIT an output format compatible with the 
COMTAJL input format has been added, allowing the representation of 
map output of any batch program on the COMTAL display unit, 

3,4* 2 Development and Maintenance Utility and Data Processing 
Programs 

The highly parametric clustering program, CLUSW, has been modified 
so that on each re clustering cycle, the thresholds on clustering distances 
and on spatial homogeneity can be changed* The statistical program 
STATS, has been provided with two additional features, the most impor- 
tant of which is the possibility of printing out the response values ox 
each point in the given test areas in all channels. The algorithm 'which 
determines the points belonging to a set of polygons describing a test 
area boundary, has been redesigned and now works correctly. 

The maximum likelihood classifier, PARAM, has been extended to al- 
low the thresholding of each class* Assignment of a point to a class will 
be inhibited if the membership probability is below the threshold. Also 
-the control cards of the PARAM program have been simplified and made 
compatible witli the other batch programs. 

Minor corrections had to be made to several batch programs. 

Since tapes of the SKY DAB S 102 13 channel- scanner became available 
to the JRC and flights with a Bendix 11 -channel scanner have been car- 
ried out for t he A GRESTE Project (see 3,4), a tape reformatting pro- 
gram, COPTRAN, has been written to convert SKY LAB tapes into full 
and half resolution and Canadian format DAJMDSAT tapes and Bendix 
tapes to the format compatible with our batch programs* Optionally 
COPTRAN performs a correction procedure for the atmospheric effects 
depending on the scan angle, through a long track averaging and smooth- 
ing technique. This procedure improved overall classification and com- 
parison performance especially for short v/ave channels, but has a ten- 
dency towards erroneous behaviour due to local inhomogeneities espe- 
cially for high contrast infrared channels* An attempt will be made to 
improve its performance by the use of low‘ order polynomial approxi- 
mation and/or reflectance gradients of homogeneous areas. 

The following features have been added to the program MAPEDIT, for 
linear geometric correction and map editing. 

The input might be any map output of the programs currently in use at 
the JRC, or any image data in line -interleaved or pixel -interleaved for- 
mat, One output format has been added (cf, 3,4*2), The coordinates of 
points in the reference map can be calculated from the corresponding 
input coordinates. 



3. 4. 3 Computer-aided Interpretation of Airborne Scanner Data 
3. 4, 3*1 MSS Flight Characteristics 

A flight was made on August 7, 1975 to acquire ct?ita from a Bendix 
M S scanner at a date on which the various rice vai'ieties would be 
in the flowering stage, considered to be a favourable phe no logical 
period for discriminating varieties of rice. 

The flight was performed two altitudes (5000 feet and 13,000 feet) 
between 9. H and 10, 3Z o'clock in the morning. The technical data of 
the seamier are reported in Table 5. 


Table 5 : BENDIX M^S scanner technical data 


Scan angle; 


O 

o 

0 


Roll-compensation: 


± 10° 


Geometric resolution 


2.5. 10'^ radian 


Band chars cterisucs: 
(Wavelengths in ^m) 

Channel Nr. 

Center 

Width 

1 

0.410 

0.06 



2 

0.465 

0.05 


3 

0.515 

0.05 


4 

0.560 

0.04 


5 

0.600 

0.04 


6 

0.640 

0.04 


7 

0.680 

0.04 


S 

0.720 

0.04 


9 

0.815 

0.09 


10 

1.015 

0.09 


11 

11.0 

6.0 


The resolution of the data is then 3. 8 ixi and 9-9 m for the altitudes 
of 5000 feet and 13, 000 feet, respectively. 

Four strips of data were acquired at both altitudes (see Fig. 7). Three 
of them are located in the AGRFSTE test site no. 1; strips 1 and Z are 
devoted to the study of rice varieties and strip 3 will be used to investi- 
gate the discrimination of poplar varieties and classes of age. The four 
strips are situated in the Po -valley region* to the west of Milan between 
the Po- and Ticino rivers. 
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Table 6 : Signification of symbols in Fig.8 and colours in Figs.S and 10 



Class 

Symbol 
in Fig.8 

Colour 

in Figs.9 & 10 

Class number 
in Table 7 


Balilla Grana Grossa 

B.G.G. 

Dark blue 

5 

CO 

03 

VP 

Gritna 

G 

Green 

6 

o 

1- 

Carnaroll 

C 

Red 

7 

0> 

u 

Rocca 

Rc 

Light blue 

8 

ir 

Arborio 

A 

Violet 

13, 14, 15 


Romeo 

Ro 

Orange 

19, 20 


Rice in general 


White 

1, 9/11, 12 


Corn 

M 

Yellow 

16, 17, IS 


Others 

///// 

Biack/dark brown 

2, 3, 4, 10 


3*4* 3,2 Discrimination of Rice Varieties!^) 

Only some results are presented here. The work will continue during 
the coming months. The data used are contained in the strip no, Z 
(see 3,.4, 3, 1), to the south of Mortara, at an altitude of 5000 feet. 

This first study covers only a' part of the prepared ground truth, which 
is sketched in Fig, 8, The classifier used -assumes euclidian distance 
between the pixels to be classed and the mean vectors of the classes of 
interest, the distance to each class being divided by the standard devia- 
tion of the class distribution. Figs, 9 and 10 show the results in a quali- 
tative way, each colour being attributed to a given class. They were ob~ 
tained by photographing the classified maps displayed on the COMTAD 
screen. The significance of colours in Figs, 9 and 10 and symbols in 
Fig, 8 is given in Table 6 and the statistical parameters of the various 
classes are collected in Table 7, 

As can be seen in Table 7 the calculations were made with 20 classes b\ 
for practical reasons some of the classes have been grouped under the 
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Table 7 : Statistical parameters of classes involved in classification process 



Clwv\. 











CiftSS 

dtiofiption 

T 


. 3 

4 

S 


? 


5 

40 

AA 

^ 67 

A 

A* 

V 

66 1 

73.770 

9,01^> 

^4.770 

19,014 

63.743 

9.386 

65.544 

14.523 

51.655 
ID. 236 

179.455 

63.043 

127,684 

33.756 

49.3QI 

98,840 

98.198 
6 . c> 47 

O.&i 

^xccMTiiiOA* 




04.679 

70.715 

56.604 

62.611 

47.416 

1 76.259 

^9,717 

32.264 

95.883 

0.63 

po\&Va fi 

jL 

V 

2* ?71 

3.693 

6.130 

6.040 

4* i'ijH 

4,201 

58,403 

38.714 

50.435 

55. OM 

4 

3 

y 


61. A2* 

69.301 

SI. 151 

66-690 

53.771 

111,177 

37,^39 

27.009 

105,26? 

€>•?? 


V 

4. I6l 

2.649 

5.004 

3,454 

4*906 

79.439 

9-LB2 

4.966 

155.948 

7.040 

sc\wifile.2 


A* 

60. 603 

C6.70L 

99,671 

85.174 

102.550 

69.929 

107.452 

08,557 

41,047 

150.636 

o.sr 

Vick I'^eUs 


V 

S.O'^l 

13.2o8 

17.022 

19.115 

5C.04S 

46.198 

60.381 

19.147 

67.672 

43.667 

S 

/■^ 

47.^07 


73,533 

58,993 

65.354 

51.089 

134,995 

05.940 

29.8B5 

104.006 

o.^g 


V 

lt>5 

H.665 

10.277 

8.510 

1^.405 

19.201 

44.180 

33.630 

61.396 

21.766 

£ 

A' 

4^. 93^ 

67.^29 

79.192 

61.5^1 

67.287 

53.548 

155.296 

114.913 

42.502 

103.172 

0.59 

flrfc- <• 6 


V 

3. 1:^3 

2. 221 

4.3D7 

13,677 

4.162 

3.4 54 

95.669 

14.735 

119-432 

6.22S 

7 

A* 

47. 44B 

6b. L9 p 

77.439 

6 1 ,b46 

63.413 

53.419 

138. 2Q6 

90,445 

27.390 

103.123 

0.€0 

- C 

V 

0. 6U 

1.260 

2.654 

18.443 

6.460 

3.079 

38.201 

3 3.240 

55.070 

0.997 



p 

r' 

8)3 

70.641 

79.777 

63.667 

71,654 

56,605 

148,341 

104.461 

35.670 

107,309 

AI6 

^\*ce-Rc 

Q 

V 

22. 691 

51.791 

67-311 

69.176 

141.414 

132.023 

IOC. 942 

64.32'* 

73.423 

76.992 


A 

V 

4 8.092 

67 ,332 

Tt,224 

60.840 

67,307 

’ 52.849 

146-233 

103. 7C9 

34.663 

104.886 

A AO 



13, 493 

.28.317 

41.2D1 

35.929 

62,725 

59.755 

156,154 

90.387 

00.157 

37.039 

Covhf^s'ife 1 

AO 


Pi, 16 5 

71.477 

81.(02 

67.492 

TS.748 

64.797 

139.856 

94,144 

33.320 

120.878 

/1.91 


V 

39, GO? 

115,331 

178.483 

130*223 

215,796 

240-951 

742*397 

169.833 

169.222 

319.070 

Hiavt 'ftce. 

41 


44, fc33 

62,181 

71.026 

57.181 

61.702 

47.375 

149.066 

106, 312 

36.245 

101,640 

o.-?2 

^iCE- 

V 

2.C2P 

Ao6l3 , 

e.463 

11.198 

1<>.629 

16.368 

28.219 

46-069 

103,457 

8-590 

BciUUcs 

Ai 

h 

44. 44ji 

65.924 

77,101 

60,751 

66.048 

51.388 

153,645 

102,947 

33.387 

103.747 

e.?! 

R»ce.- 

V 

1. 733 

2.790 

6.961 

4.259 

10.657 

6.644 

37.736 

22,374 

34.800 

6-346 

j^owiq 


V 

50,743 

0,492 

71.675 

0,C3l 

B3.728' 

1.447 

65.860 

1.629 

71.953 

2.12L 

57.510 

2.010 

166.010 

5,532 

111.546 

6.920 

39.203 

49.243 

107.065 

6.231 

O.Sf 

fiv'ce, A 



4 3, 676 

02.692 

70,17 3 

56,629 

61 .4 4 0 

46.7 24 

150.965 

112.786 

38.942 

103.332 

0,44 

^tce. A 

4A' 

3. 236 

1.573 

3.97 L 

4.O02 

2.476 

3.115 

6.106 

8,694 

47.236 

9.490 

t 


:C 

59,717 

70.2J5 

G3.4L3 

65.828 

7-1.79 2 

57.113 

161.526 

117.231 

44.901 

no. 366 

0.S6 

t?ice- A 

AS 

0. 539 

1,410 

3.453 

2.460 

2.557 

2.219 

18,548 

12.192 

49,219 

4.951 


46 

/ 

V 

5 0, 357 
C.$39 

69,698 

2.306 

81,364 

6.016 

64,743 

3.930 

71.591p 

5.757 

57.916 

7.763 

177.469 

137,377 

118.590 

S5.341 

46.574 
ICO. 305 

1L3.163 

8.S87 

O.s? 

scavvxpls-l 

4? 


5 4.607 

73.971 

84.250 

69.^32 

81.715 

66,487 

145.274 

97.549 

35*243 

120.769 

o:n 

Co'fW 

4. 294 

l2,lo9 

ia, 0 L 6 

21.274 

45,406 

33,836 

68.331 

38,00? 

66,639 

21,355 

S6tvW]ple^ 

4? 


52. 634 

66.500 

73.524 

61.050 

71.327 

55.639 

121.347 

92.068 

31.764 

114.556 

0,42 


C. R60 

l.llO 

1.55? 

20.643 

2,620 

2.827 

4.162 

2.616 

30.776 

10.543 


49 


44, 139 

*62.984 

72.602 

56.345 

59.093 

44,389 

153.031 

109.522 

39-039 

95,003 

o.gr 


V 

C.793 

1.023 

1,624 

1.096 

1.143 

1.414 

9.672 

12.45B 

659,897 

6.057 


^0 


41.9B3 

1.300 

62.323 

3.7Q5 

73.133 

12.271 

59.‘t99 

6.634 

65.330 

6.983 

52.106 

6.453 

147. 7DC 
136.629 

97.271 

43.640 

29.036 

42.457 

103.H90 

6.836 

0.64 

Bi'cc 


U ; mean value of each channel over the proper training sample 

V : variance of the same data 

2 ~ ; weighted standard deviation for each class, where a = 

1=2 Ml 

* Rica varieties BGG, G, C, Rc end A sampled together 
** M^er fields, com and uncultivated areas sampled together 
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same coloar on Figs, 9 and 10, In paxticular the white colour covers 
all the rice varieties which were not sampled on the area processed 
here (i. e. classes no. 1, 11 and 12 in Table 7) and b composite rice 
class - no. 9 - composed by sampling together the rice varieties sam- 
pled in the area, i. e. BGG, G, C, Rc and A), The black/dark brown 
colour, on the other hand, covers a composite class made up by sam- 
pling together wet meadows, corn and uncultivated zones (class no. 10), 

and poplars and wet meadows present as distinct classes (class no. 2. 3 
and 4). ' 

The sampling was made on small portions of the fields labelled with the • 
appropriate symbols on the area processed (Fig. 8} or outside the area 
for classes no. 1, 11 and 12, 

The data from chaniiel 1 were not used, due to an unacceptable level of 
noise. The calculations were made without any membership threshold 
for the classes, therefoi-e no point was left unclassified. 

In Fig. 9 the data were processed without correcting the channel values 
for the atmospheric effect due to the scan angle (100° angle from edge to 
edge of the swath). In Fig. 10 the data were corrected from this effect 
with the so-called '‘long track averaging” technique with, as a result, 
a better discrimination between rice and corn on the one hand, and rice 
and ^^others^' on the other hand* 

An overall examination of the classifications displayed show a partial 
discrimination of rice varieties with somewhat better results in Fig. 10. 
Variety BGG is rather well-recognized but outside its limits is confilsed* 
with variety C; partial confusion occurs between varieties G and A where- 
as Rois poorly recognized and Rc is not recognized at all. 

Corn discrimination is also rather good in Fig. 10 (with atmospheric cor- 
rection). 

The above comments on rice varieties discrimination are obviously valid 
only for those rice fields whose variety is identified on the ground truth. 
The other rice fields, not labelled on the ground truth, contain unidenti- 
fied varieties or a mixture of varieties. It must also be noted that 

concision occurs between rice and non- rice, outside the rice region to 
the right of tlie a rea studied. 

These qualitative results must be considered to be preJiminary. It is hoped 
to make improvements by further refining the atmospheric correction tech- 
nique cur renay under study, the sampling procedure and the classification 
algorithm. A quantization of the ground truth will enable us on the other 
hand, to present quantitative results under the form of a confusion matrix. 
On the other hand, the data from the higher flight (13, 000 feet) will also 
be used and the effect of the loss in spatial resolution will be investigated. 




Fig. 3 : Ground-truth for rice varieties 

Aoprox. scale Vi5.000 



Fig. 9 : Classification results for rice varieties without 

atmospheric correction on the data 






Fig. 8 : Ground truth for rice virieties 

Aop*'ox. scale 5.000 


I 



Fig. 10 : Classification results for rice 

varieties after correction of the 
data from atmospheric effect 


,1 
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P 
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3.4,4 Computer-aided Interpretation of Satellite Data, Classification 
of Natural Beech Forests 

Two LANDSAT scenes were free of clouds over the studied region, 
namely those of July 3, 1975 and September 13, 197 5. 

These two scenes were treated first by level-slicing on channel 7 in 
order to facilitate the recognition of geographic details. 

Supervised methods of classification have been tlie most widely used. 
Signatures have been determined by making statistics over homogeneous 
regions. 

Two classed of beeches were considered. One in sunlight and one in 
shadow. Even though chestnuts were of no direct interest to the AGRESTE 
Project, they had to be included in the claSbification because their sig- 
nature is so close to that of beeches that misclassification could easily 
occur betw^een these two species. Meadows covering the bottom of the 
valleys were also included. Automatic classification was then produced 
by different methods based either on Euclidean distance: a pixel is clas- 
sified in the **nearest” class, or on the maximum likelihood principle: 
the pixel is classified in the class to which it has the highest probability 
to belong. 

As the ground truth was not exhaustive (more than 50% of the surface is 
not characterized), the classification methods must use thresholds i, e. 
limits of distance or probability over or under which a pixel is not clas- 
sified. 

The best results were obtained by a maximum likelihood method (PARAM 
program): the results from July 3 are displayed in Fig, 12 and those of 
September 1 3 in Fig. 13, 

Qualitative comparison with the geographic map reproduced schematic- 
ally in Fig, 11 i^ quite good, and most of the beech forests are correctly 
identified. 

In the scene of July 3, it is possible to separate beeches (red) from chest- 
nuts (blue). Small isolated chestnut woods in the valley-bottom meadows 
are also identified. 

The scene of September 13 confirms the results but here the discrimina- 
tion between beeches and chestnuts is not so good, probably due to a less 
favourable period of vegetation. Fig. 14 shosvs a comparison of the two 
scene -classifications. The comparison is made on entirely deciduous 
forests i, e. beeches + cheotnuts, blue corresponds to the classification 
of July3, red to the classification of September 13. Purple, the super- 
imposition of these two colours, gives the intersection of the two classi- 
fications i. e. the pixels classified in the same way on both scenes. It 
can be seen that correspondence is quite good. 


0 > 


I 


i 

! 

i 

i 


I 


i 

i 

i 

I 



Fig. 11 : Ground truth reported from aerial photos 

red: beeches 
blue: chesnuts 
green: meadows 
white: not characterised 

The river "Stura di Demonte" and the villages of Demonte (east) 
and Aisone (west) are reported to facilitate localisation 
Scale Is about 1 cm = 1.2 km 
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Fig. 12 : 


Results of parametric classification 

red: beeches 

blue: chesnuts 

green: meadows 

black: not classified 

The frames on the figure correspond 

to those of ground truth 

Scale is about 1 cm = 2.2 km 
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Fig. 13 : Results of parametric classification of 
September 13th 

red: beeches 
blue: chesnuts 
green: meadows 
black: not classified 



Fig. 14 : Comparison of classifications of June 3rd 
and September 13th 
Comparison is made on total deciduous forest : 
purple: zone classified on both scenes 
blue: zone classified only on June 3rd 
red: zone classified only on September 13th 
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Due to tlie lack of an exhaustive ground truth, quantitative compari- 
sons are not easy. In fact, the level of the threshold used for each 
class is one of the important parameters for the final population of 
the class. Up to now, the thresholds have been fu-:ed empirically taking 
into account the variance-covariance matrix of each class. 
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